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Abstract. In 1976, Dodziuk and Patodi employed Whitney forms to define 
a combinatorial codifferential operator on cochains, and they raised the ques- 
tion whether it is consistent in the sense that for a smooth enough differential 
form the combinatorial codifferential of the associated cochain converges to 
the exterior codifferential of the form as the triangulation is refined. In 1991, 
Smits proved this to be the case for the combinatorial codifferential applied 
to 1-forms in two dimensions under the additional assumption that the ini- 
tial triangulation is refined in a completely regular fashion, by dividing each 
triangle into four similar triangles. In this paper wc extend Smits's result to 
arbitrary dimensions, showing that the combinatorial codifferential on 1-forms 
is consistent if the triangulations axe uniform or piecewise uniform in a cer- 
tain precise sense. We also show that this restriction on the triangulations is 
needed, giving a counterexample in which a different regular refinement pro- 
cedure, namely Whitney's standard subdivision, is used. Further, wc show 
by numerical example that for 2-forms in three dimensions, the combinatorial 
codifferential is not consistent even for the most regular subdivision process. 



1. Introduction 

Let M be a n-dimensional polytope in M" , triangulated by a simplicial complex 
Th which we orient by fixing an order for the vertices. (Although wc restrict to 
polytopes for simplicity, several of the results below can be easily extended to 
triangulated Riemannian manifolds.) Wc denote by A'= = h^{M) the space of 
smooth differential fc- forms on M. The Eiiclidean inner prodiict restricted to M 
detc!rminc!s the Hodge star operator A*^ A"^*^, and the inner product on A*^ given 
by ((i, (,') = J u A The space L^A*^ is the completion of A*^ with respect to this 
norm, i.e., the space of differential /c- forms with coefficients in . We then define 
HK*' to be the space of forms u in L/^K^ whose exterior derivative du, which may be 
understood in the sense of distributions, belongs to L^A'^^^. These spaces combine 
to form the LF' de Rham complex 

-> i?A° -4 iJ A^ A • • • -4 if A" ^ 0. 
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Viewing the exterior derivative d as an unbounded operator L^A'^ to L^A'^+^ with 
domain HA'' , we may define its adjoint d* . Thus a differential fc-form u belongs to 
the domain of d* if the operator v {u, dv)i^2j^k is bounded on L^A'^"^, and then 

{d*u,v) i^2/^k-i — {u,dv)i^2j^k, V e HA''^^. 

In particular, every u which is smooth and supported in the interior of Af belongs 
to the domain of d* and d*u = -kd-ku. 

Let Ak{Th) denote the set of /c-dimensional simplices of Th- We denote by Ck{Th) 
the space of formal linear combinations of elements of I^k{Th) with real coefficients, 
the space of fc-chains, and by C''{Th) — Ck{Th)* the space of fc-cochains. The 
coboundary maps d'^ : C'^{Th) — > C''^^{Th) then determine the cochain complex. 
The de Rham map Rh maps A*"' onto C''{Th) taking a differential fc-form u to the 
cochain 

(1.1) Rhu : Ck{Th) ~> R, 

The canonical basis for C''{Th) consists of the cochains Ot, t € Ak{Th), where Ur 
takes the value 1 on r and zero on the other elements of Afc(7^). The associated 
Whitney form is given by 

k 

Wha-r = kl ^{-iy\i dAo A . . . A (iAi A . . • A d\k, 

i=0 

where Aq, . . . , are the barycentric coordinate functions associated to the vertices 
of the simplex listed in order (thus Xi is the continuous piecewise linear function 
equal to 1 at the ?th vertex of r and vanishing at all the other vertices of the 
triangulation) . The span of Whar, t e Ak{Th), defines the space of A^ of Whitney 
fc-forms. Its elements are piecewise affine differential fc-forms which belong to HA'' 
and satisfy dAjj C Af^'^^ . Thus the Whitney forms comprise a finite dimensional 
subcomplex of the de Rham complex called the Whitney complex: 

0^A° ^A^^... Aa^'^O. 
The Whitney map Wh maps C''{Th) isomorphically onto Ajj and satisfies 

(1.2) Whd'c^dWhC, ceC''{Th), 

i.e., is a cochain isomorphism of the cochain complex onto the Whitney complex. 
Although Whitney fc-forms need not be continuous, each has a well-defined trace 



on the simplices in Ak{Th), so the de Rham map (1.1) is defined for u e A^. 
The Whitney map is a one-sided inverse of the de Rham map: RhWhC = c for 
c G C''{Th)- The reverse composition tt/i — WhRh '■ A'' A^j defines the canonical 
projection into A^. 

In |3j and [4], Dodziuk and Patodi defined an inner product on cochains by 
declaring the Whitney map to be an isometry: 

(1.3) {a,b) = {Wha,Whb)L2Ak, a,beC''iTh). 

They then used this inner product to define the adjoint S'' of the coboundary: 

(1.4) {6^a, b) = (a, #6), a, b e C'^iTh)- 

Since the coboundary operator d'^ may be viewed as a combinatorial version of the 
differential operator of the de Rham complex, its adjoint S'^ may be viewed as a 
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combinatorial codifferential, and together they define the combinatorial Laplacian 
on cochains given by 

The work of Dodziuk and Patodi concerned the relation between the eigenvalues of 
this combinatorial Laplacian and those of the Hodge Laplacian. 

Dodziuk and Patodi asked whether the combinatorial codifferential S"^ is a con- 
sistent approximation of d* in the sense that if we have a sequence of triangulations 
Th with maximum simplex diameter tending to zero and satisfying some regularity 
restrictions, then 

(L5) lim WWhS^RhU ~ d*u\\ ^ 0, 

h 

for sufficiently smooth u e A*^ belonging to the domain of d* . Here and henceforth 
the norm |j • || denotes the norm. 

Since C^{Th) and A^ are isometric, we may state this question in terms of Whit- 
ney forms, without invoking cochains. Define the Whitney codifferential : A| — 
by 

(L6) (d*u,w)i2Afe-i = (u,dw)i2A<c, u e hi, V e A|;-i. 

Combining ([L2]), ([O]), and ([O]), we see that d* = WhS^Wf^^. Therefore, WhS^Rh = 
d'^TTh, and the question of consistency becomes whether 

(1.7) lim lld^vT/iU — d*u|| = 0, 

h 

for smooth u in the domain of d* . 



In Appendix H of PJ, the authors suggest a counterexample to (1.71 for 1-forms 
(i.e., A: = 1) on a two-dimensional manifold, but, as pointed out by Smits [7j, the 
example is not valid, and the question has remained open. Smits himself consid- 
ered the question, remaining in the specific case of 1-forms on a two-dimensional 
manifold, and restricting himself to a sequence of triangulations obtained by regular 
standard subdivision, meaning that the triangulation is refined by dividing each tri- 
angle into four similar triangles by connecting the midpoints of the edges, resulting 
in a piecewise uniform sequence of triangluations. See Figure [S] for an example. In 
this case, Smits proved that (1.5) or, equivalently, (1.7) holds. 

Smits's result leaves open various questions. Does the consistency of the 1-form 
codifferential on regular meshes in two dimensions extend to 

• Mesh sequences which are not obtained by regular standard subdivision? 

• More than two dimensions? 

• The combinatorial codifferential on /c- forms with fc > 1? 

In this paper we show that the answer to the second question is affirmative, but the 
answers to the first and third are negative. More precisely, in Section [2] we present 
a simple counterexample to consistency for a quadratic 1-form on the sequence of 
triangulations shown in Figure [T] While these meshes are not obtained by regu- 
lar standard subdivision, they may be obtained by another systematic subdivision 
process, standard subdivision as defined by Whitney in [51 Appendix II, § 4]. Next, 
in Section [Sj we recall a definition of uniform triangulations in n-dimensions which 
was formulated in the study of superconvergence of finite element methods, and we 
use the superconvergence theory to extend Smits's result on the consistency of the 
combinatorial codifferential on 1-forms to n-dimensions, for triangulations that are 
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uniform or piecewise uniform. In Section|4] we provide computational confirmation 
of these results, both positive and negative. Finally, in Section [5] we explore nu- 
merically the case of 2-forms in three dimensions and find that the combinatorial 
codifFerential is inconsistent even for completely uniform mesh sequences. 

2. A COUNTEREXAMPLE TO CONSISTENCY 

We take as our domain M the square (—1,1) x (—1,1) C M^, and as initial 
triangulation the division into four triangles obtained via drawing the two diagonals. 
We refine a triangulation by subdividing each triangle into four using standard 
subdivision. In this way we obtain the sequence of crisscross triangulations shown 
in Figure [ij with the mth triangulation consisting of 4™ isoceles right triangles. 
We index the triangulation by the diameter of its elements, so we denote the mth 
triangulation by Th where h = 4/2*". Using this triangulation, the authors of [5] 
showed that superconvergence does not hold for piecewise linear Lagrange elements. 



Figure 1. 72, 71, 71/2, 7i/4, the first four crisscross triangulations. 

Define p : M ^ M by p(a;, y) ^ x - x^/S and let u = dp = (1 - x^)dx £ A^{M). 
Now for q G HhP{M) (i.e., the Sobolev space H^{M)), we have 

f dq , dq , \ dq , dq , 
-kdq = * -^dx + —dy = —dy - T^-dx, 
\ox ay J ox ay 

so 

{u,dq) L2/^i — f uA-kdq= f {1 — x^)-^ dx dy — f 2xqdxdy — {2x,q)]^2j^o. 
Jm Jm ox 7 m 

Thus u belongs to the domain of d* and d*u = 2x. As an alternative verification, 
we may identify 1-forms and vector fields. Then u corresponds to the vector field 
(1 — x^, 0) which has vanishing normal component on dM, so belongs to the domain 
of d* = - div and d*u = - div(l - x'^,0) = 2x. 

Set Wh = dl^TTfiU. Now Wh G A'l, i.e., it is a continuous piecewise linear function. 
The projections tt/i into the Whitney forms form a cochain map, so tt/jM — iThdp — 
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d-TThP — grad iThP where iThP is piecewise hnear interpolant of p. Thus Wh G A° is 
determined by the equations 



(2.1) 



Whq dx dy 



M 



/ grad TThp ■ grad qdxdy, g G . 
Jm 



It turns out that we can give the solution to this problem explicitly. Since Wh is a 
continuous piecewise linear function, it is determined by its values at the vertices of 
the triangulation Th- The coordinates of the vertices are integer multiples of h/2. 
In fact the value of Wh at a vertex (x, y) depends only on x and for /i < 1 is given 

by 



Wh{x,y) = < 



— /i, X = —1, 

0, — l<a;<l, xa multiple of h, 

h, X = 1, 

-6 + 2/1, x = -l + h/2, 

6x, —1 + h/2 < X < 1 — h/2, X an odd multiple of h/2, 

6-2/1, x = l-h/2. 



A plot of the piecewise linear function Wh is shown in Figure [3]for /i = 1/2. To 
verify the formula it sufhces to check (2.1) for all piecewise linear functions q that 



vanish on all vertices except one. There are several cases depending on how close 
the vertex is to the boundary, and the computation is tedious, but elementary. Here 
we only give the details when the vertex is {x,y) with — l + ft,/2 < x < 1 — h/2 
and x is an odd multiple of h/2. To this end, let q be the piecewise linear function 
that is one on vertex (x, y) and vanishes on all the remaining vertices. In this case, 
the support of q is the union of the four triangles Ti,T2,T^,T4^ that have {x,y) as 
a vertex (see Figure [5]) . According to the formula, in the support of q, one has 
Wh = Qxq. A simple calculation then shows that the left hand side of (2.1 1 is 



/ Whqdxdy = 6xy / q^dxdy — Axm, 



where m — /i^/4 = jTij for any i. 



To calculate the right-hand side 


of 


(|2.1[) for this q, we 






'(1,0), 


on Ti , 


grad q = 




(0,1), 


on Ti, 




(-1,0), 


on T3, 






1(0,-1), 


on T4. 



and 



grad 7r;,p = 



•(p(x)-p(x-|),0), 
(lb(x+|)-p(x-|)],p(x) 
(p(x+|)-p(x),0), 

. (1 b(x + I) - p(x - I)] , 1 b(x + I) + p{x 



;[p(a 



-p(^-t)]), 
|)]-P(^)), 



onTi, 
on T2, 
on T-i, 
on Ti^. 



Hence, 
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/ grad TThP ■ grad qdxdy — / T^hP • grad q dx 



dy 



16 
ft2 



— ) + pyx + —)\)m — 4:xm. 



This verifies (2.11 for this piecewise hncar function q. 





Finally, we note that, since Wh essentially oscillates between 6a; and 0, it does 
not converge in to d*u (or to anything else) as h tends to zero. 
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3. Consistency for 1-forms on piecewise uniform meshes 

We continue to consider a sequence of triangulations Th indexed by a positive 
parameter h tending to 0. We take h to be equivalent to the maximal simplex 
diameter 

ch < max AiavaT<Ch, 

for some positive constants C, c independent of h (throughout we denote by C and c 
generic constants, not necessary the same in different occurences). We also assume 
that the sequence of triangulations is shape regular in the sense that there exists 
c> 

piT) > cdiamT, 

for all T E Th and all h, where p{T) is the diameter of the ball inscribed in T. 

We begin with some estimates for the approximation of a fc-form by an clement 
of A^. For this we need to introduce the spaces of differential forms with coefficients 
in a Sobolev space. Let to be a non-negative integer and u a /c-form defined on a 
domain M C M", which we may expand as 

(3.1) u= u,^...,^dx'' A ■ ■ ■ Adx'>'. 

Using multi-index notation for partial derivatives of the coefficients Ui-^...i^. , we define 
the mth Sobolev norm and seminorm by 

l<ji <---<ifc <n |Q|<m 

|u|^„^fc = X/ X/ \\^°''^ii--ik\\l^{M)i 

l<!l<---<'ifc<n |a|=m 

and define the space iJ™A'^(M) to consist of all fc-forms in M for which the Sobolev 
norm ||u||ifmAfc is finite. 

With these notations, we can state the basic approximation result that for any 
shape regular sequence of triangulations there is a constant C such that 

(3.2) inf < C/i||u||^iA'., uEH^h^{M). 

For a proof, see [U Theorem 5.8]. Since H^A^ is dense in L^K^ , this implies that 

(3.3) dist(/,A^) := inf II/- w|| ^ as /i^ 0, f e L^A''{M). 



In addition to the best approximation estimate (3.2), we also need an 0{h) 
estimate on the projection error ||u — 7r/iu||. For this we require more regularity of 
u, since tt^m is defined in terms of traces of u on fc-dimensional faces, which need 
not be defined on H^A^. 

Lemma 3.1. Let {Th\ be a shape regular sequence of triangulations of M C M" 
and k an integer between and n. Let £ be the smallest integer so that £ > (n — 
k)/2. Then there exists a constant C , depending only on n and the shape regularity 
constant, such that 



(3.4) \\7ThU-u\\L2Ak <C Y.h'^HH^A', ueH'A\M). 

m— 1 
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Proof. First we note that the canonical projection is defined simplex by simplex, 
as 

{TThU)\T = TTt{u\t), 

where, for v a /c-form on T, tttv is its interpolant into the space of Whitney forms 
on the single simplex T. Therefore, it is enough to prove that 

e 

(3.5) \\u - tttuWl^^.^t) <CY. /i">lff"AMT), n e H'A''{T), 

m— 1 

with the constant C depending on T only through its shape constant. We prove 
this first for the unit right simplex in M", T, with vertices at the origin and the n 
points (1,0, ... ,0), (0, 1,0, . . .), .... Since £ > {n — k)/2, we obtain, by the Sobolev 
embedding theorem, that |!7i'j.M||^2Afc(f ) — ^ll'"llHfA'»(f )' ^^"^ triangle 
inequality, 

Now let u = nl J^u, a constant fc-form on T equal to the average of u. Then 

TTfU = U, SO 

||u-7r^u||^2Afc(f) = II(m-u) -7r^(u-u)||^2Afe(f) 

i 

< C\\u — u\\HiA>'{f) — ^(11" ~ ^11 L2A'=(f) + ^ l^l//"iA'=(f ))' 

rn — 1 

where we have used the fact that u is a constant form, so its mth Sobolev seminorm 
vanishes for m > 1. Now we invoke Poincare's inequality 

Putting things together, and writing u instead of u, we have shown that 

e 

(3.6) ||u-7r^u||^2A'=(f) ^ C'Xl l^lff'"A'=(f)' ueH^A''{f). 



This is the desired result (3.51 in the case T — T. 

To obtain the result for a general simplex, we scale via an afSne diffeomorphism 
F -.f ^T. If u is the fc-form on T given by ([sl]), then 

(3.7) F*u^ E i^n-^.oF)^---^dx^^A---Adx^''. 

{l<ii<...<ifc<n} ji,... jfc=l 

Each of the partial derivatives / dx^'' is a constant bounded by h. Using the 
chain rule and change of variables in the integration, we find that 

(3.8) c|F*u|^„^.(^) < (voir)-i/2 I^^I^^AHT) < C|F*u|^„A.(^), 
where the constants c and C depend only on m and n and the shape regularity 



constant of T. Combining iS.Gl and (3.8) we get 



||m - 7rTu||_f^2Afc(T) < C(voir)^/^/i - 7r^u||^2Afc(f) 



< C(voir)l/2/j-'= \^\H^A>^if) <CJ2 h'^'lulH^A^T), 

m— 1 m— 1 



which estabhshes (3.5 1. □ 
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Our approach to bounding the norm of the consistency error is to relate it to 
another quantity which has been studied in the finite element literature, namely 

„^ ... - TThU, dVh) 

(3.9) Ah{u) := sup — . 



Theorem 3.2. Assume the approximation property (3.3). Then, for any smooth 
u G L^A'^ belonging to the domain of d* we have 

lim — d^7r/iu|| limA;i(u)=0. 

h h 

This follows immediately from the following lemma. 

Lemma 3.3. Let 1 < k < n, and let u G L?K^ he smooth and in the domain of d* . 
Then 

(3.10) Ah{u) < \\d*u - dl-KhuW < dist(d*u, A^"^) + Ah{u). 
Proof. The first inequality is straightforward. For any Vh G A^~^ , 

hh\\ hh\\ 
For the second inequality, we introduce the L^-orthogonal projection Ph : L'^K*'^^ — > 
A^~^ and invoke the triangle inequality to get 

(3.11) \\d*u^dl-Khu\\ < \\d*u-Phd*u\\ + \\Phd*u-dlTrhu\\ = dist(d*u, A^"^) + ||w|i, 
where w — Phd*u — d*f^'KhU G A^. Now 

(3.12) WwW"^ = {Phd*u - dlTThU,w) = {u- tthu, dw), 
hence 

3.13 W = r-- < sup r r = Ah{u), 

which completes the proof. □ 

Thus we wish to bound {u — i^hU^dyh) /\\vh\\ for smooth u in the domain of d* 
and Vh G A^. An obvious approach is to apply the Cauchy-Schwarz inequality and 
then use the approximation estimate (3.4) 

(3.14) \{u-TThU,dvh)\<\\u-TThu\\\\dvh\\ <Ch\\u\\HtKk\\dvh\\. 

To continue, we need to bound ||c?Wft||/||uh|| for Vh an arbitrary nonzero element 
of A^. Because Ajj consists of piecewise polynomials, it is possible to bound its 
derivative in terms of its value using a Bernstein type inequality or inverse estimate. 
This gives that 

(3.15) \\dvh\\<Ch-^\\vhl Vh^Al 

where h = minxgA„(T/i) diamT. Unfortunately, even if we assume that our trian- 
gulations are quasiuniform, i.e., that h>ch for some fixed c > 0, this just leads to 
the bound 

Ah{u) < C||u||^«Afc, 

which does not tend to zero with h. In fact, we cannot hope to get a bound which 
tends to zero without further hypotheses, since, as we have seen, even for the nice 
mesh sequence and form u considered in the previous section, d'^ is not consistent, 
and so Ah{u) does not tend to zero. 
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Nonetheless, for very special mesh sequences it is possible to improve the bound 



(3.141 from first to second order in h. This was established by Brandts and Kffzek 
in their work on gradient superconvergence [5]. The mesh condition is embodied 
by the following concept. 

Definition 3.4 ([2]). A triangulation T on M is called uniform if there exist n 
linearly independent vectors ei, . . . , e„, such that 

(1) Every simplex in T contains an edge parallel to each ej. 

(2) // an edge e is parallel to one of the Cj and is not contained in dM , then 
the union of simplices containing e is invariant under reflection through 
the midpoint of e, i.e., 2me ~ x £ Pe for all x G Pe- 

The crisscross triangulations shown in Figure [T] satisfy the first condition of the 
definition, but not the second, and so are not uniform. On the other hand, the mesh 
sequence obtained starting from a single triangle, or from a division of a square into 
two triangles and applying regular standard subdivision, is uniform. See the first 
two rows of Figure |4] A uniform triangulation of the cube in n dimensions is 
obtained by subdividing it into to" subcubes, and dividing each of these into n\ 
simplices sharing a common diagonal, with all the diagonals of the subcubes chosen 
to be parallel. The 3-D case is shown in Figure |4) We refer to for more details. 




Figure 4. Uniform triangulations. 

Theorem 3.4 of [5] claims that if {Th} is a shape regular family of uniform 
triangulations of M, and if u is a smooth 1-form, then there exists a constant 
C > such that 

\{TThU - u,dvh)\ < Ch'^\\u\\H-^A^\\dvh\\, 

for all Vh € A-ln H'^(M) and h > 0. Here H'^{M) denotes the space of H^(M) 
functions with vanishing trace on dM. However, their proof uses the inequality (cf. 
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(1.5) of [2 ) 

\\T^hU - w||l2ai < C /l||u||_H-iAi, 

where C is a constant independent of u. This would imply that tth can be continu- 
ously extended to i/^A^, which is impossible for n > 3. Fortunately, the proof in [2] 
works verbatim, if the above inequality is replaced by (3.4). Hence, the following 
result is essentially proved in [5]. 

Theorem 3.5. Let {Th\ be a shape regular family of uniform triangulations of M , 
and let u be a smooth 1-form. Furthermore, let £ be the smallest integer so that 
£ > (n — l)/2. Then there exists a constant C > such that 

(3.16) \{7:hu - u,dvh)\ < Ch^\\u\\H.^i\\dvhl 
for all Vh£ Aln H^{M) andl> h>0. 

Next we consider piecewise uniform sequences of triangulations. 

Definition 3.6. A family Th of triangulations of the polytope M is called piecewise 
uniform if there is a triangulation T of M such that for each h, Th is a refinement 
ofT and for each T € T and each h, the restriction ofTh to T Cz A„(T) is uniform. 

If, as in [7j , we start with an arbitrary triangulation of a polygon and refine it by 
standard regular subdivision, the resulting sequence of triangulations is piecewise 
uniform. This is illustrated in Figure [5] The following theorem shows that d* 
is consistent for 1-forms on piecewise uniform meshes, thus generalizing the main 
result of [7] from 2 to n dimensions. 

Theorem 3.7. Assume that the family of triangulations {Th} is a shape regular, 
quasiuniform, and piecewise uniform. Let u S H^A^{M) be a 1-form in the domain 
of d* , where £ is the smallest integer satisfying £ > (n ~ l)/2. Then we have 

(3.17) lim d;,7r,,u|| = 0. 

Proof. Let T denote the triangulation of M with respect to which the triangulations 
Th are uniform. We will apply Theorem|3.5|to the uniform mesh sequences obtained 
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by restricting Th to each T To this end, let K — UteT denote the skeleton 
of T, and set 

^h = \J{TeTh\TnK^(d}. 

We can decompose an arbitrary function Vf,. G as 

(3.18) Vh = Wh+Y,vl, 

where Wh G A° is supported in 'Sh and G A° is supported in T. Indeed, we just 
take Wh to coincide with v at the vertices of the triangulation contained in K and 
to vanish at the other vertices, while v]^ = v at the vertices in the interior of T 
and vanishes at the other vertices. Because the mesh family is shape regular and 
quasiuniform, there exist positive constants C, c such that 

c\\vr<h J2 Hx)\'<C\\vr, veAl 

from which we obtain the stabihty bound 

(3.19) \\wh\\ + Y.\\^^\\^^\\^'^1 

T<£T 



Using the decomposition (3.18) of Vh we get 

I {iThU - u, dvh) I < I (tt/jU - u, dwh) I + X! I ~ '^'"h) I 

Ter 

< Ch\\u\\HtA_^(Y.H)\\dwh\\ +Ch^ ^ \MH'^h^(M)\\dvl\ 
(3.20) Ter 

< C\\u\\H'-K^(^^)\\wh\\L^-{T.^) + ChY ||m||h^ai(M)I|w/T|I 

Ter 

< C (||M||/f«Ai(Sfc) + /i||"llHfAi(M)) hhW-, 

where we have used the Cauchy-Schwarz inequality, the projection error estimate 
(3.4), the second order estimate (3.16) (which holds on the uniform meshes on each 
T), the inverse estimate of (3.15), and the L^-stability bound (3.19). Since the 
volume of E/j goes to as /i 0, so does ||w||_f/«Ai(Sfc)- Thus Ah{u) vanishes with 



/i, and the desired result is a consequence of Theorem 3.2 □ 



Remark 3.8. The preceding proof shows that as long as the triangulation is mostly 
uniform, in the sense that the volume of the defective region goes to as /i — > 0, we 
obtain consistency. One can also extract information on the convergence rate. For 
instance, using the fact that is 0{h), we obtain ||it||/f«Ai(s,,) ^ C\/h\\u\\Qejyi for 
wGC^Ai(M). 



4. Computational experiments for 1-forms 

In this section, we present numerical computations confirming the consistency 
of for 1-forms on uniform and piecewise uniform meshes in 2 and 3 dimensions, 
and other computations confirming its inconsistency on more general meshes. The 
four tables of the section display the results of computations with various mesh 
sequences. In each case we show the maximal simplex diameter h, the number of 
simplices in the mesh, the consistency error ||(i^7r/i/ — and the apparent order 
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inferred from the ratio of consecutive errors. All computations were performed 
using the FEniCS finite element software library [5]. 

The first two tables concern the problem on the square described in Section [2] 
i.e., the approximation of d*u where u = (1 — x^)dx. Table [l] shows the results 
when the piecewise uniform mesh sequence shown in the Figure [5] is used for the 
discretization. Notice that the consistency error clearly tends to zero as 0{h). 



n 


h 




triangles 


error 


order 


1 


S.OOe- 


-01 


20 


6.25e- 


-01 




2 


2.50e- 


-01 


80 


3.08e- 


-01 


1.02 


3 


1.25e- 


-01 


320 


1.56e- 


-01 


0.98 


4 


6.25e- 


-02 


1,280 


7.85e- 


-02 


0.99 


5 


3.12e- 


-02 


5,120 


3.94e- 


-02 


1.00 


6 


1.56e- 


-02 


20,480 


1.97e- 


-02 


1.00 



Table 1. When computed using the 2-dimensional piecewise uni- 
form mesh sequence of Figure [5] the consistency error tends to 
0. 



By contrast. Table [2] shows the counterexample described analytically in Sec- 
tion [2) using the mesh sequence of Figure [T] obtained by standard subdivision. 
In this case, the consistency error does not converge to zero, as is clear from the 
computations. 



Table 2. 
error does not tend to 0. 



n 


h 




triangles 


error 


order 


1 


5.00e- 


-01 


16 


1.15 




2 


2.50e- 


-01 


64 


1.50 


-0.38 


3 


1.25e- 


-01 


256 


1.60 


-0.09 


4 


6.25e- 


-02 


1,024 


1.62 


-0.02 


5 


3.12e- 


-02 


4,096 


1.63 


-0.01 


6 


1.56e- 


-02 


16,384 


1.63 


-0.00 



Similar results hold in 3 dimensions. We computed the error in d^it on the 
cube (—1, 1)'^ where again u is given by (1 — x^)dx. We calculated with two mesh 
sequences, both starting from a partition of the cube into six congruent tetrahedra 
all sharing a common edge along the diagonal from (— 1,— 1,— 1) to (1,1,1). We 
constructed the first mesh sequence by regular subdivision, yielding the meshes 
shown in Figure [6| These are uniform meshes, and the numerical results given 
in Table [3] clearly demonstrate consistency. For the second mesh sequence we 
applied standard subdivision, obtaining the sequence of structured but non-uniform 
triangulations shown in Figure [7] In this case d*^ is inconsistent. See Table [4] 

5. Inconsistency for 2-forms in 3 dimensions 

We have seen that for, 1-forms, d\ is consistent if computed using piecewise 
uniform mesh sequences, but not with general mesh sequences. It is also easy 
to see that consistency holds for n-forms in n-dimensions for any mesh sequence. 
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Figure 7. As in 2D, the mesh sequence in 3D obtained by stan- 
dard subdivision is not uniform. 



This is because the canonical projection tt/i onto the Whitney 71-forms (which are 
just the piecewise constant forms) is the orthogonal projection. Now if Vh 
is a Whitney (n — l)-form, then dvh is a Whitney n-form, so the inner product 
{u — TThU^dvfi) = 0. Thus Afi{u), defined in (3.91, vanishes identically, and so 
is consistent by Theorem |3.2[ Having understood the situation for 1-forms and n- 
forms, this leaves open the question of whether consistency holds for A:-forms with 
k strictly between 1 and n. In this section we study 2-forms in 3 dimensions and 
give numerical results indicating that is not consistent, even for uniform meshes. 

Let M = (1 - x^){l - y^)dx A dy, a 2-form on the cube M = (-1,1)^. The 
corresponding vector field is (0, 0, (1 — x^){l — y^)) which has vanishing tangential 



ON THE CONSISTENCY OF THE COMBINATORIAL CODIFFERENTIAL 15 



n 


h 


tetrahedra 


error 


order 


1 


l.OOe+00 


48 


1.69e+00 




2 


5.00e-01 


384 


9.70e-01 


0.80 


3 


2.50e-01 


3,072 


5.13e-01 


0.92 


4 


1.25e-01 


24,576 


2.63e-01 


0.96 


5 


6.25e-02 


196,608 


1.33e-01 


0.98 


6 


3.12e-02 


1,572,864 


6.69e-02 


0.99 



Table 3. The consistency error for d^^ on 1-forms in 3-D tends to 
zero, when using the uniform mesh sequence of Figure [6j 



n 


h 


tetrahedra 


error 


order 





l.OOe+00 


48 


l.Slc+OO 




1 


5.00e-01 


384 


2.71e+00 


-0.58 


2 


2.50e-01 


3,072 


3.02e+00 


-0.16 


3 


1.25e-01 


24,576 


3.11e+00 


-0.04 


4 


6.25e-02 


196,608 


3.13e+00 


-0.01 



Table 4. The consistency error for on 1-forms in 3-D using the 
nonuniform mesh sequence of Figure [7] does not tend to zero. 



components on dM . Therefore u belongs to the domain of d* and d*u is the 1- 
form corresponding to curlit, i.e., d*u = —2(1 — x^)ydx + 2x{l — y^)dy. Table [5] 
shows the consistency error ||fijj7r^u — d*u\\i,-2{t^i computed using the sequences of 
uniform meshes displayed in Figure |6] This mesh sequences yields a consistent 
approximation of d* h for 1-forms, but the experiments clearly indicate that this is 
not so for 2-forms. 



n 


h 




triangles 


error 


order 


1 


1.00c4 


-00 


48 


1.59e-f00 




2 


5.00e- 


-01 


384 


1.18e-f00 


0.43 


3 


2.50e- 


-01 


3072 


l.OOe-^OO 


0.24 


4 


1.25e- 


-01 


24576 


9.47e-01 


0.08 


5 


6.25e- 


-02 


196608 


3.37e-f00 


-1.83 



Table 5 . The consistency error does not tend to zero for 2-forms, 
even on a uniform mesh sequence. 
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